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Abstract. The finite Heisenberg group knows when the dimension of Hilbert space is a square number. Remarkably, it then 
admits a representation such that the entire Clifford group — the automorphism group of the Heisenberg group — is represented 
by monomial phase-permutation matrices. This has a beneficial influence on the amount of calculation that must be done to 
find Symmetric Informationally Complete POVMs. I make some comments on the equations obeyed by the absolute values of 
t-O 1 the components of the SIC vectors, and on the fact that the representation partly suggests a preferred tensor product structure. 
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THE HEISENBERG GROUP 

Q-L The name "Heisenberg group" for the group we are interested in became widely used some 50 years after Heisenberg 

did the path breaking work that Arkady Plotnitsky describes elsewhere in this volume — and which is indeed closely 
related to this group. The very foundations of quantum mechanics are closely related to the Heisenberg group. 

Actually there are many Heisenberg groups. They all admit a defining representation in terms of upper triangular 
matrices of the form 




(1) 



> 
ON 

Such matrices form a nilpotent group as soon as the entries x, p, and <p can be added and multiplied together in such 
a way that they form a group under addition. In other words, the entries must belong to a ring, but the ring can be 
chosen freely. If the entries are real numbers we obtain a three dimensional Lie group, whose Lie algebra includes 
the commutator between position and momentum — figuring so prominently in Heisenberg's discoveries. If the entries 
belong to the ring of integers modulo N we obtain a finite Heisenberg group, denoted H(N), which has an essentially 
unique unitary irreducible representation in a Hilbert space of N dimensions. Other options are available, such as 
choosing the entries to belong to a finite field. This leads to finite groups represented in //-dimensional Hilbert spaces, 
where p is a prime number. Unless N is a prime number the latter groups differ from H(N). 

I have no intention to describe this zoo of groups in detail here. Suffice it to say that it played a major role in two 
of the most brilliant episodes of human thinking, nineteenth century geometry yl and twentieth century physics — 
and that I believe they will play a major role in this century too. The focus in my talk was on the role they play in the 
problem of finding Symmetric Informationally Complete POVMs, or SICs for short. In the laboratory SICs correspond 
to a special kind of (doable [2]) measurements, of interest for the foundations of quantum mechanics JU. My personal 
motivation for working on the problem will be revealed at the end. 

A SIC is a collection of A^ 2 unit vectors in a Hilbert space of dimension N, such that they form a resolution of the 
identity, i.e. a POVM, and such that 

|< W | W )| 2 =-^(1+JV-S„), 1<I,J<N 2 . (2) 

Finding such collections of vectors is surprisingly hard. Zauner's conjecture states that they do exist for any N, that they 
can be chosen such that they form an orbit of the Heisenberg group H(N), and such that every vector in the collection 
is left invariant under a special order three element of the automorphism group of H(N) Jj]. This automorphism group 
is what we call the Clifford group here. The evidence for Zauner's conjecture comes from explicit constructions, not 
from any structural understanding. Scott and Grassl give a wonderful summary of the status of this conjecture in the 
concluding section of their recent paper . 



It is known that H(N) is the only group that can serve the purpose if N is a prime not equal to 3 la]. It is also known 
that the analogous problem concerning equiangular lines in real Hilbert spaces does not have a solution in general I2D. 

There is a problem with a similar flavour, which is that of finding a complete set of Mutually Unbiased orthonormal 
Bases, or MUB. Here it is believed that a solution exists only if the dimension equals a prime number or a power of 
a prime number. The existing solutions [8] are again closely related to the Heisenberg groups, but then one uses the 
version based on finite fields, which is not the same as H(N) unless the dimension N is a prime number. It is odd that 
the SIC and MUB existence problems differ in this way. They also differ in that — at least from one point of view — the 
latter is partly understood: it is known that the existence of a complete set of MUB is equivalent to the existence of 
a unitary operator basis forming a flower with N + I petals, in the sense that the A^ 2 operators in the basis can be 
divided into N+l petals of mutually commuting operators, such that two different petals have only the unit element in 
common |9|]. If we restrict ourselves to unitary operator bases whose elements form a group, this means that we must 
look for groups with a sufficient number of maximal abelian subgroups. It is known that no group can do better than 
the Heisenberg groups in this regard lloll . (Why one should restrict oneself to unitary operator bases of group type is 
not clear at all. But this is beside the point for the moment.) 



SQUARE DIMENSIONS 

It would seem as if H(N), the Heisenberg group over the ring of integers modulo N, has a kind of "universal" structure 
that is insensitive to the choice of N, except for some special features that emerge when N is prime. However, closer 
inspection reveals that there is something very special about square dimensions too. 
To see this, recall that the usual presentation of H(N) is in terms of the root of unity 

Zni 

co = eT . (3) 

and generators X and Z subject to the relations 

ZX = coXZ X N = Z N = 1 . (4) 

When presented like this the Heisenberg group admits a unique unitary representation — up to unitary equivalence [11]. 
The usual choice then is to use the representation where the cyclic subgroup generated by Z is diagonal. But suppose 
that the dimension N = n 2 isa square number. Then it is readily seen that 



Z n X n =X n Z n . (5) 



This gives rise to a preferred maximal abelian subgroup all of whose N elements are of order n. What we [12] refer 
to as the phase permutation basis is the basis in which this special subgroup is diagonal. Its basis vectors are labelled 
\r,s), where r and s are integers modulo n, and we define the phase factor 

q = e " = CO . (6) 
With the phase conventions we used the representation of the group generators is then 

\r,s+l) ifi+l^Omodn 
Xlr.s} = { ' (7) 
q r \r,0) ifi+l=0modn 

Z\r,s) =CO s \r-l,s) . (8) 

This leads to 



X"\r,s)=q r \r,s) , Z"\r,s) = q s \r,s) . (9) 

All matrices here are monomial unitaries, also known as phase-permutation matrices, meaning that exactly one element 
in each column is non-zero, and equal to a phase factor, and similarly for the rows. Actually this much is true for the 
usual representation too. The reason why we refer to the new basis as the phase-permutation basis lfl2ll is that its use 
has in its train that the entire Clifford group is given by phase-permutation matrices — and this is not at all true for the 
usual representation of the Heisenberg group. Recall that the Clifford group consists of all unitary operators U such 



that UX'ZW^ is again an element of the Heisenberg group — and that the Clifford group plays a very important role 
when one tries to understand the SIC problem, which is our main concern for the moment. 

To show that the Clifford group is represented by phase-permutation matrices it is enough to observe that, being 
an automorphism group, it must permute the maximal abelian subgroups — known as stabiliser groups in quantum 
information theory — among themselves. But it must also preserve the order of all group elements. When the dimension 
equals n 2 there exists a unique stabiliser group consisting solely of elements of order n. This is the stabiliser group that 
determines the new basis. The Clifford group can only permute its elements among themselves, and the remarkable 
property of the basis follows immediately. 

There is an interesting connection to the theory of theta functions here; in particular the 'new' representation is used 
by Mumford in his presentation of the theory of theta functions from the point of view of the Heisenberg group 111 311 . 
The Heisenberg group over the real numbers can be represented by generators S x and T y acting on the space of analytic 
functions through 

(S x f)(z) = f(z + x) , (T y f)(z) = e^+^fiz + zy) , (10) 

where z is a complex number with positive imaginary part. If we restrict ourselves to the lattice of integer valued x 
and y there is a unique function invariant under these transformations, namely the theta function 

*=-oo 

Using a net with somewhat finer meshes we catch an invariant n 2 -dimensional vector space of functions. This vector 
space admits a natural basis consisting of theta functions of rational characteristic. Thus we define a subgroup of the 
real Heisenberg group by insisting that na and nb are integers. Let 0o,o(z) — 9(z). As a basis for the vector space we 
take the functions 



e il ,b(z)=S b T a 9(z) , na,nbeZ. (12) 
The phase-permutation representation is recovered from the observations that 



S c 9 a ,b = 9 a ,b+c ) T c e ab = e 27tibc 9 a+cb , na,nb,nceZ (13) 



and 



a +x,b+y = e 2%iay e a ,b , x,y,na,nb € Z . (14) 

The reader may be familiar with the case n = 1/2, discussed in Whittaker and Watson |3; in a notation closer to 
Jacobi's we have (0oo> 9 n i , © i n , 0i i ) = (O3, 64, 02, 9i). Via elliptic curves the theta functions play a prominent role 

u >2 2 ,u 2 '2 

in the geometry of complex projective space, so the connection may prove helpful for the SIC problem, or in some 
other way — although at the moment this is a vague suggestion only. 

In another direction, note the temptation to write \r,s) = \r) ® \s). Eq. (O then says that the special subgroup we are 
basing the representation on is represented by local operators, that is by operators of the form U\®U2- Moreover the 
generator Z is a local operator — as indeed it was in the representation proposed by Weyl ifTHl . where Z is diagonal. If 
then the basis vectors are separable one finds 

Z = diag(Lo),...,a) Ar - 1 ) =diag(l ) $,...,$")®diag(l,fl>,. ...ffl" -1 ) . (15) 

The unitary that connects the two representations is a local operator with respect to this tensor product structure, so this 
property is preserved. We can conclude, for instance, that the n 4 vectors in a given SIC must fall in n multiplets with 
the same degree of entanglement in each. (Curiously, even more is true for some N = A SICs: all the 16 vectors have the 
same concurrence 11511 . I do not understand why.) However, note that our two requirements — that \r,s) forms a phase- 
permutation basis and that \r,s) = \r) <E) \s) — contains some ambiguities. The Main Property of the phase-permutation 
basis is that it leads to a representation of the automorphism group of the Heisenberg group by phase-permutation 
matrices, and this property is kept if the basis is changed by an arbitrary unitary phase-permutation matrix. In general 
this will disrupt the tensor product structure. Still the phase-permutation basis partly determines the tensor product 
structure of the Hilbert space, and this may turn out to be interesting — but again this is a somewhat vague idea at the 
moment. 



THE SIC EQUATIONS FOR THE ABSOLUTE VALUES 



I advertised that the phase-permutation basis simplifies the SIC problem to some extent. The main reason is that it is 
well adapted to exploit Zauner's conjecture, according to which any SIC vector is left invariant by a special element of 
the Clifford group. This does help, although not as much as one could hope for. The details can be found in the paper 
that we wrote ll 211 — and since the paper appeared very recently I cannot add much to the story here. 

What I can add are a few comments that were left out of our paper. In a previous Vaxjti meeting Marcus Appleby 
discussed a subset of the SIC equations ||2), which concern only the absolute values of the components of the SIC 
vectors 11611 . He assumed that the SIC is an orbit of the Heisenberg group. In this case the SIC is completely determined 
by a single fiducial vector | y/o) sucn that 

r i i=j=o 

|(V/ |rz^|v/o)| 2 = < (16) 
I nTT otherwise . 

The equations involve quartic polynomials in N variables, and most known solutions have been derived using Grobner 
bases and dedicated computer programs such as MAGMA [5]. Still it is interesting to manipulate them a bit further, in 
order to see what is involved. Let the components of the fiducial vector be given by 

Za = sfFJ^ • (17) 

In the usual representation it is known that the SIC equations are equivalent to the equations ifTvL PTsIl 



JV— l 9 

*>« = *TT' (18) 



N-1 



a=0 



£ = ^y, xe{l,2,...N-l}, (19) 



a=0 



N-1 



Y, ZaZa+k-iZu+kZa-i = , Z, k ^ . (20) 

a=0 

Interestingly there are N equations that do not involve the phases at all. Moreover these equations have an interesting 
geometrical interpretation, arising when we project all vectors onto the simplex spanned by the basis vectors in a cross 
section of the convex body of density matrices. The equations imply that the image of the A^ 2 SIC vectors is a regular 
simplex of a characteristic size, inscribed in the larger simplex 11611 . Unfortunately the geometric interpretation also 
suggests that the absolute values pi cannot be determined by these equations alone, since the orientation of the image 
simplex is not determined. 

If we examine eqs. (118H19t we see that they they remain the same if the sign of the integer x is changed. Hence there 
are at most k + 1 equations independent equations, regardless of whether N = 2k or N = 2k + 1 . We also see that they 
imply 

/N-1 \ 2 

(21) 

So normalisation is included — as in fact it was from the start, before the SIC equations were brought to this form. If 
N = 2k we can also derive that 

/, . . . \ 2 

' (22) 





N+l 



which again is easily imposed, and means that the equations can be linearised if N — 2. However, I cannot find any 
further relation that strikes the eye and is easy to handle. 

Very similar equations can be deduced in the phase permutation basis, where the absolute values are naturally 
denoted by y/p7 s - Recall that the dimension is N = n 2 . From the equations for a Heisenberg covariant SIC one deduces 
that fTT 



n-J 

Y, PrsPrs 

r,s=0 

n-1 

Y PrsPr+x.s 
r,s=0 



N+l 



1 



N+l 



(23) 



(24) 



Here x,y are integers modulo n, not both zero. In itself this is a slight improvement on the previous result; the number 
of independent equations goes up to k + 2 if N = 2k. This means that the absolute values are determined by these 
equations alone if N = 4. 

Let N = 4. Then there are four equations for the absolute values, all of which can be brought to a form analogous to 
eqs. d2Tl22b : 

I 
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(25) 



After taking the square roots we obtain a linear system of equations. Since the components appear symmetrically in 
the equations we may assume without loss of generality that poo > poi > pio > p\\ > 0. This determines three of the 
signs. If the sign of the fourth square root is positive one finds that 



Poi =Pio = Pii = 



5 - V5 
20 



(26) 



The negative sign is inconsistent with p rs > 0. Hence the moduli are completely determined, and solving for the 
phases afterwards is not too hard. Notice that the calculation is not only straightforward, it also gives every Heisenberg 
covariant SIC in four dimensions, with no input from numerical searches. 

The phase-permutation basis comes into its own when we impose Zauner's conjecture. By some further permutations 
and rephasings of the basis vectors one can block diagonalise a Zauner unitary, so that it consists of three by three 
blocks of the form 




(27) 



together with some diagonal elements. Each such block contributes three eigenvalues ( 1 , e 2m / 3 ( e 4 ^'/3^ ( j-jjg spectrum. 
Given that the spectrum of a Zauner unitary is known |01 one can deduce that there will be 3 diagonal elements if 
= 3/, and only one if N = 31 + 1. Regardless of whether the block diagonalisation is carried through, one concludes 
that an invariant vector will have / + 1 independent components only, where N = n 2 = 31 or 31 + 1 . The resulting 
equations (123H24I ) for the absolute values are then easy to manage, at least for = 9, 16,25 which we did explicitly. 
But it turns out that the number of independent equations is too small to determine all the absolute values — and solving 
the equations that involve the phases is not easy at all lfl2ll . In fact it has never been done, in any basis, for N — 1 6 , 25 Q 



CONCLUSION 



The phase-permutation basis partly suggests a preferred tensor product structure Jff" ® Jff" of the Hilbert space 
Jf? N . The commuting operators X" and Z" are local operators. This is interesting in itself. It also has the wonderful 



This is what I said in my talk. Eventually things worked out better than I expected — the final version of our paper does contain the solution for 
.V I60|. 



consequence that the entire Clifford group is represented by monomial unitary matrices. This simplifies the SIC 
existence problem, but more ideas are needed to trivialise it in dimensions N > 4. Perhaps some can be found. 

All of this reminds me of Fermat's Last Theorem. The statement of Fermat's theorem is very simple. In itself it 
appears to be a curiousity only, but in the course of proving it mathematicians were led to develop deep theories that 
changed the face of their subject. I like to think that the SIC existence problem has a similar status within quantum 
mechanics. 
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